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Abstract 

In  ionic  clusters  the  polarization  depends  on  the  dynamics  of  charge 
delocalization.  The  polarization  follows  the  charge  motion  when  the 
characteristic  time  of  motion  is  much  longer  than  the  polarization  time 
(a^^batic  polarization).  For  the  opposite  case  (diabatic  polarization),  the 
polarization  is  determined  by  an  average  field.  These  two  kinds  of 
polarization  are  considered  for  the  ionic  clusters  Xe*R  and  benzene+R 


(R-He,Ne,Ar) . 
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1.  Introduction 

Ionic  atoms  and  molecules  or  neutral  molecules  in  excited  ionic  states 
can  form  clusters  by  attracting  neutral  closed- shell  atoms  or  molecules.  If 
the  ionization  potential  of  the  attracted  particles  is  high  enough,  they  will 
remain  neutral  inside  the  cluster.  As  examples  of  such  clusters  we  can 
mention  Ar^Ar^  [1-4],  aromatic  cations  with  rare  gas  atoms  [5-7],  Cl  Xej,  or 
ClXe^Xe^in  a  rare  gas  matrix  [8,9].  These  clusters  are  kept  together  mainly 
by  polarization  forces. 

In  simple  model  calculations  the  energy  of  an  ionic  cluster  is 
presented  as  a  sum  of  repulsive,  dispersive  and  polarization  components. 
Whereas  the  repulsive  and  dispersive  components  can  be  considered  in  the 
pairwise  approach,  the  polarization  energy  depends  significantly  on  the 
cluster  geometry  and  the  charge  distribution.  The  last  point  is  of  importance 
for  the  case  of  a  molecular  ionic  core  of  a  cluster  where  the  charge  is 
delocalized  between  some  number  of  atoms.  The  charge  delocalization  in  a 
positive  ionic  cluster  core  (molecule)  can  be  described  as  a  motion  of  a 
positive  hole  with  some  characteristic  time  of  the  hole  "jump"  between 
adjacent  atoms.  This  time  is  of  the  order  of  W/E,  where  E  is  the  coupling 
energy  responsible  for  charge  delocalization.  Since  neutral  atoms  and/or 
molecules  also  need  some  time  to  be  polarized,  the  interactions  between  them 
and  the  delocalized  charge  is  influenced  by  dynamical  effects.  When  the 
polarization  time  is  much  smaller  than  the  characteristic  time  of  the  hole 
motion  (adiabatic  polarization),  then  the  polarization  will  follow  the  hole 
motion  and  the  polarization  energy  will  become  an  average  of  the  polarization 
energies  of  all  hole  locations.  Thus  in  the  case  of  one  neutral  atom 
polarized  by  an  ionic  molecule,  the  adiabatic  polarization  in  the  point  dipole 
approximation  is 
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U  - 

a 


(1) 
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where  a  is  Che  neutral  atom  polarizability,  a^  is  the  hole  population  at  the 
i-th  atom,  and  F ^  is  the  field  generated  by  a  hole  located  at  i-th  atom.  For 
the  opposite  case,  where  the  polarization  time  is  much  larger  than  the 
characteristic  time  of  hole  motion  (diabatic  polarization) ,  the  neutral  atom 
will  be  polarized  by  an  average  field, 

"a  -  •  !  •  <2> 

i 

As  an  example  of  a  system  (unrealistic,  unfortunately)  which  demonstrates  a 

dramatic  difference  between  two  kinds  of  polarizations,  we  can  mention  an 

ionic  ring  with  a  neutral  atom  at  the  ring  center.  For  the  adiabatic  case  the 

2 

polarization  energy  is  -af\/2  (any  i)  ,  whereas  for  the  diabatic  case  it  is 
zero. 

The  aim  of  this  work  is  to  present  the  polarization  energy  for  the 
general  case  of  arbitrary  relation  between  the  polarization  time  and  the 
charactistic  time  of  the  hole  motion.  The  dynamical  effects  of  the  hole 
motion  on  the  polarization  energy  will  be  described  below  in  the  terms  of  the 
correlation  between  the  hole  location  and  the  polarized-atora  electrons. 

2.  Polarization  Energy 

Let  us  consider  a  cluster  consisting  of  a  positively-charged  molecular 
ion  M  and  a  neutral  atom  A  polarized  by  the  molecular  ion.  The  positively- 
charged  hole  of  the  isolated  molecular  ion  is  described  by  the  wave  function 


u 


where  is  the  atomic  orbital  of  the  hole  located  at  the  i-th  atom  of  the 
molecular  ion  M.  In  the  zero  overlap  of  atomic  orbitals  (ZOAO)  approximation, 


<^i^j>  “  5ij  •  }  ai  “  1  • 


the  energy  of  the  molecular  ion  M  is 

40)-h-2^  }  aiajhij  ' 

i  i  j>i 


where  is  the  energy  of  diabatic  state  with  located  at  the  i-th  atom  hole 
and  h^j  is  the  interatomic  exchange  energy.  For  the  sake  of  simplicity  the 
neutral  atom  is  assumed  to  be  of  S-symmetry,  and  the  wave  function  of  the 
polarized  atom  is  described  as  the  superposition  of  S-  and  P-states: 


*r<v 


bsS<rA>  +  }  bkPk<cA>  ■ 


k  -  x,y,z 


bs  *  K 


It  is  assumed  that  the  polarization  is  weak,  Ib^j  «  bg. 
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We  shall  present  the  wave  function  of  the  M-A  cluster  in  a  form  which 
takes  into  account  the  correlation  between  the  hole  location  and  the  neutral 
atom  state: 


♦<¥«> 


‘Cs,iS(rA> 


Ck.lW> 


k-1 


(7) 


We  have  assumed  that  the  polarization  is  weak  so  that  the  charge  (hole) 
distribution  in  M  is  not  affected  by  the  polarized  atom  A: 
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Taking  into  account  these  relations  in  the  wave  function  (7),  one  obtains 


*(rMrA> 
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In  the  ZOAO  approximation  this  wave  function  provides  the  energy  expression 
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(10) 


where  E^  ^  and  E^  are  the  energies  of  the  separated  molecular  ion  M  and 
neutral  atom  A,  respectively,  Is  the  energy  of  interaction  between  the 
molecular  ion  and  non-polarized  atom  A,  and  U  is  the  polarization  component  of 
the  M-A  interaction, 
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In  Eq.  (11)  €  is  the  atomic  S-»P  excitation  energy,  and  ^  is  the  k  -  x,y,z 
component  of  the  interaction  between  the  located  at  the  i-th  atom  charge  and 
the  polarized  atom.  In  the  point  dipole  approximation  [Eqs.  (l)-(2)]  the  ion 
and  the  atom  polarized  by  it  are  supposed  to  not  overlap  with  one  another. 
The  neglect  of  overlap  leads  to  the  unphysical  divergence  of  ^  at  small 
separations.  In  order  to  prevent  this  divergency,  a  damping  function  7  is 
introduced  into  the  polarization  energy  expression  (4,7,10], 


Vi  •  hi  Fk.i  “ 


where  ^  is  the  k-x,y,z  component  of  the  field  generated  at  the  center  of 
the  atom  A  by  the  hole  charge  located  at  the  i-th  atom  of  the  molecular  ion  M, 
and  n  is  the  transition  dipole  moment. 


a  -  e<P  Sx>  -  e<P  yS>  -  e<P  zS>. 
x  y  z 


(13) 


In  the  two-state  atomic  model  applied  here  (Eq.  (6)],  the  transition  dipole 
moment  can  be  expressed  by  the  atom  polarizability  a  as 


(13') 
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The  damping  function  7  is  taken  from  Ref.  4  as 


■*1  +  M12r1/3 


t<«)  -  (1  +  [^]  ] 


where  r  is  the  separation  between  the  ion  i  and  the  atom  A,  and  p ^  and  p  are 
the  atomic  radii  of  the  ion  and  atom,  respectively. 

Applying  the  variational  principle  to  the  polarization  energy  (11),  we 
obtain  for  the  coefficients  ^  the  system  of  linear  equations 


t€  +  a  ^  aihij1Ck,i  '  5  ’  *^.1  * 


Let  us  consider  the  solution  of  Eqs.  (15)  for  two  extreme  cases,  e  » 
|h^  |  and  <  «  |h.. |.  For  the  first  case  the  coefficents  ^  are 


C.  .  -  W,  .  . 

k.i  €  k.i 


Substituting  these  coefficients  into  Eq.  (11)  and  taking  into  account  Eqs. 
(12)  and  (13'),  one  obtains 


V  -  Ua  -  -f-  }  4  n  }  4,1  -  ■  2}  al1'l<?l)2  • 


which  is  the  expression  for  the  adiabatic  polarization  energy.  It  differs 
from  expression  (1)  by  the  damping  functions  7^  only.  Physically  it  is  clear 
that  the  condition  e  »  |h^  |  leads  to  the  adiabatic  polarization,  since  it 
means  that  the  polarization  time  is  much  smaller  than  the  characteristic  time 
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of  the  hole  motion.  For  the  opposite  case  of  diabatic  polarization  «  « 
|hij  |  ,  the  polarization  time  is  large  so  that  the  neutral  atom  "feels"  the 
average  field.  The  solution  of  Eqs.  (15)  in  this  case  is 


a  ,  -  ^  ^  a2  ».  , 

k,i  «  L  J  kJ 
j 

V  -  ud  -  -I  <5  *Vh  V 


According  to  this  expression,  the  average  field  in  the  diabatic  polarization 
is  weighted  by  square  roots  of  the  damping  functions  7^. 

For  the  particular  case  of  a  homonuclear  diatomic  ion  (i-1,2),  the 
solution  of  Eqs.  (15)  is 


Ck,i  “ 


(t+h)wk_1+huk_2 


(c+h)Wk  ^hWk  l 


Jl  £(€+2h)  ’  k’2  Jl  £(£+2h) 


which  provides  the  polarization  energy 


u  -  •  !'<vi  *  V2> +  ?i +  ^2  ?2)2' 


The  adiabatic  (£  »  h)  and  diabatic  (e  «  h)  polarization  energies  are 


Ua  *  '“■Vl  +  72F21  ' 


Ud  *  -S^  +  Jh  ?2>: 


The  polarization  energy  (21)  becomes  equal  to  the  average  of  Ua  and  when 
h  —  0.5e.  It  follows  that  for  the  case  of  a  diatomic  ion  the  polarization  is 
adiabatic  for  h  «  0.5e  and  diabatic  for  h  »  0.5e. 

Dealing  with  ionic  clusters,  we  need  to  know  whether  to  use  the 
adiabatic  (17)  or  diabatic  (19)  expression.  For  most  of  ionic  clusters,  the 
adiabatic  condition  «  »  h  is  fulfilled  since  the  exchange  term  h  is  usually 
of  the  order  of  1-2  eV,  whereas  the  excitation  energy  e  of  rare  gas  atoms,  for 
example,  is  of  the  order  of  10  eV.  It  is  of  interest  also  to  know  how  large 
the  differences  are  in  the  adiabatic  and  diabatic  polarization  energies.  In 
order  to  study  these  differences,  we  shall  consider  below  the  clusters  formed 
by  Xe+  ion  and  rare  gas  atoms  and  by  an  aromatic  cation  and  rare  gas  atoms. 

3.  Xe*R  (R-He.Ne.Ar)  Clusters 

In  a  simple  model  calculation,  the  Xe*-R  potential  is  presented  as  a 
sum  of  the  pairwise  van  der  Waals  (vdW)  interaction  and  the  R-atom 
polarization  energy.  We  shall  assume  that  the  vdW  interaction  between  the 
neutral  R  atom  and  the  charged  Xe  atom  is  the  same  as  between  two  neutral 
atoms,  Xe-R.  We  shall  describe  the  vdW  potential  by  the  L-J(12-6)  expression, 


U 


vdW 


£[(f)12 


(*)6] 


(24) 


We  have  used  empirical  values  [11-14J  for  the  energy  E  and  distance  a.  In  the 
expression  for  the  damping  function  (14)  the  sum  of  atomic  radii  is  taken 
equal  to  the  vdW  parameter  a, 
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The  Xe^R  clusters  are  of  T-shaped  geometry  (Fig.  1).  The  Xe^-R  potentials  for 
this  geometry  are 

Ud(p)  -  2UvdW(r)  -  C7(r)/r4  (25) 

Ud(p)  -  2UvdW(r)  -  C7(r)cos29/r4  (26) 

C  -  ae2/2  ,  (27) 

where  r  -  p/cosd  is  the  Xe-R  separation  and  p  is  the  distance  from  R  atom  to 
the  center  of  Xe*R  clusters.  Thus,  the  distances  of  zero  potential  p^  and 
minimum  energy  p^  are  smaller  in  the  adiabatic  than  the  diabatic  approach  by 
0.03  A  and  0.03-0.05  A,  respectively  (Table  1).  The  dissociation  energy  in 
the  adiabatic  approach  is  larger  by  some  9-13%. 

4 .  Benzene  -R  Clusters 

The  ionization  potential  of  benzene-rare  gas  clusters  is  red-shifted 
relative  to  that  of  isolated  benzene  due  to  the  polarization  of  the  rare  gas 
atom  by  ionized  benzene  [5],  In  model  calculations  [5-7]  the  rare  gas  atom 
polarization  was  considered  in  the  averaged  field  of  delocalized  benzene+ 
charge,  i.e.,  by  the  diabatic  approach.  According  to  our  estimations, 
however,  the  adiabatic  approach  has  to  be  applied  in  this  case.  Taking  into 
account  that  in  the  benzene-R  and  benzene+-R  clusters  the  rare  gas  atom  R  is 
located  on  the  benzene  symmetry  axis  Cg,  the  adiabatic  and  diabatic 
polarization  energies  differ  one  from  another  by  the  coefficient  cos  6,  like 

in  Eqs.  (25)-(26).  If  to  take,  for  example,  the  benzene+-He  cluster,  the 
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minimum  energy  distance  is  -  3.0  A  [5]  which  gives  cos  0  -  0.82  and  the 
difference  between  the  adiabatic  and  diabatic  polarization  energies  of  some 
0.004  eV.  In  the  region  of  zero  potential  (Pq  -  2.6  A,  r q  -  3.0  A,  cos  9  - 
0.78)  the  difference  between  the  two  polarization  energies  is  0.006  eV.  As  in 
the  case  of  Xe*R  clusters,  such  a  difference  in  the  polarization  energy  is  not 
crucial  for  our  results,  but  it  is  also  not  negligible,  so  that  in  further 
benzene  +  -R  clusters  calculations  it  is  desirable  to  use  the  adiabatic 
polarization  energy. 

5 .  Conclusions 

In  ionic  clusters  formed  by  a  polyatomic  ionic  core  and  a  neutral  atom, 
the  atomic  polarization  is  usually  faster  than  the  charge  motion,  so  that 
adiabatic  polarization  takes  place  and  the  polarization  energy  is  determined 
as  an  average  over  the  charge  locations.  For  the  opposite  (diabatic)  case, 
the  polarization  energy  is  determined  by  an  average  over  the  field  of  the 
charges.  As  shown  for  simple  ionic  clusters,  the  difference  between  the 
adiabatic  and  diabatic  polarization  energies  is  noticeable  but  not  large. 
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Table  1.  Comparison  of  the  parameters  of  Xe£R  clusters  (Fig.  1)  obtained  with 

the  adiabatic  and  diabatic  polarization  energies  [Eqs.  (25) - (26) ] .  p ^  is  the 

distance  of  zero  potential,  p  is  the  distance  of  minimum  energy,  and  D  is 

the  dissociation  energy.  The  distances  (o,p  ,p  )  are  in  A,  and  the  energies 

o  m 

(E,D  )  are  in  eV. 
e 

adiabatic  diabatic 


R 

a 

E 

Ref. 

'0 

^m 

D 

e 

po 

^m 

D 

e 

He 

3.56 

0.00964 

12 

3.00 

3.50 

0.0105 

3.03 

3.55 

0.0095 

Ne 

3.42 

0.0258 

13 

2.87 

3.35 

0.0138 

2.90 

3.39 

0.0122 

Ar 

3.80 

0.0654 

14 

3.27 

3.80 

0.0674 

3.30 

3.88 

0.0620 

Figure  Caption 

Arrangement  of  atoms  in  Xe+R  cluster. 


